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We present theoretical predictions for electron scattering on oxygen and calcium isotopic chains. 
The calculations are done within the framework of the distorted-wave Born approximation and 
the proton and neutron density distributions are evaluated adopting a relativistic Dirac-Hartree 
model. We present results for the elastic and quasi-elastic cross sections and for the parity-violating 
asymmetry. As a first step, the results of the models are tested in comparison with some of the 
data available for elastic and quasi-elastic scattering on ^^O and ''"Ca nuclei. Then, the evolution 
of some nuclear properties is investigated as a function of the neutron number. We also present a 
comparison with the parity-violating asymmetry parameter obtained by the PREX Collaboration 
on ^''^Pb and give a prediction for the future experiment CREX on ''^Ca. 
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I. INTRODUCTION 

The nuclear response to an external probe is a pow- 
erful tool to investigate the structure of hadron systems 
such as atomic nuclei and their constituents. In par- 
ticular, electron-scattering reactions have provided the 
most complete and detailed information on nuclear and 
nucleon structure [1-6]. Electrons interact with nuclei 
predominantly through the electromagnetic interaction, 
which is well known from quantum electrodynamic and 
is weak compared with the strength of the interaction 
between hadrons. The scattering process is therefore ad- 
equately treated assuming the validity of the Born ap- 
proximation, i.e., the one-photon exchange mechanism 
between electron and target. The virtual photon, like the 
real one, has a mean free path much larger than the tar- 
get dimensions, thus exploring the whole target volume. 
This contrasts with hadron probes, which are generally 
absorbed at the target surface. Moreover, the ability to 
vary independently the energy uj and momentum q of 
the exchanged virtual photon transferred to the nucleus 
makes it possible to map the nuclear response as a func- 
tion of its excitation energy with a spatial resolution that 
can be adjusted to the scale of processes that need to be 
studied. 

Several decades of experimental and theoretical work 
on electron scattering have provided a wealth of infor- 
mation on the properties of stable nuclei. Nuclear charge 
density distributions and charge radii have been deter- 
mined from the analysis of elastic electron scattering 
data [7, 8]. At low energy loss bound nucleons may 
produce excited states as a result of single-particle (s.p.) 
transitions and/or collective motion. Spectroscopic anal- 
ysis of such excitation mechanisms is the traditional field 
of nuclear physics. The systematic investigation of inelas- 



tic electron scattering has given the necessary support to 
the foundation of many-body theories applied to the nu- 
clear systems. Spin, parities, and the strength and struc- 
ture of the transition densities connecting the ground and 
excited states have been studied [2, 9]. In comparison 
with hadron scattering, which also gives important in- 
formation, only electron scattering can be used to de- 
termine the detailed spatial distributions of the charge 
transition densities for a variety of single-particle and 
collective transitions. At energy loss little above parti- 
cle emission threshold the quasi-bound giant resonances 
occur. Here the possibility of independently varying oj 
and <7, in conjunction with the coincident detection of 
nucleons emitted in the decay process, allows a detailed 
analysis of the various types of collective motions that are 
responsible for contributions of different multipolarities 
to such resonances. Several electric and magnetic giant 
multipole resonances have been discovered and studied 
using electron scattering [10, 11]. 

At energies above giant resonances a large broad peak 
occurs at about lo = g^/(2mAr), where uin is the nu- 
cleon mass. Its position corresponds to the elastic peak 
in electron scattering by a free nucleon. It is quite natu- 
ral to assume that a quasi-free process is responsible for 
such a peak with a nucleon emitted quasi-elastically. If 
the nucleons were indeed free, the peak would be sharp 
and would just occur at lo ~ / {2mN), corresponding 
to the energy taken by the recoiling free nucleon. A shift 
in the position of the peak is produced by the nuclear 
binding, while a broadening of the peak is produced by 
Fermi motion. 

Coincidence (e, e'p) experiments in the quasi-elastic 
(QE) region represent a very clean tool to explore the 
proton-hole states. A large amount of data for the ex- 
clusive (e, e'p) reaction have confirmed the assumption of 



2 



a direct knockout mechanism and has provided accurate 
information on the s.p. structure of stable closed-shell 
nuclei [4, 5, 12-17]. The separation energy and the mo- 
mentum distribution of the removed proton, which allows 
to determine the associated quantum numbers, have been 
obtained. From the comparison between experimental 
and theoretical cross sections it has been possible to ex- 
tract the spectroscopic factors, which give a measurement 
of the occupation of the different shells and, as a conse- 
quence, of the effects of nuclear correlations, which go 
beyond a mean field description of nuclear structure. 

In the inclusive (e, e') process only the scattered elec- 
tron is detected and the final nuclear state is undeter- 
mined, but the main contribution in the region of the 
QE peak still comes from the interaction on single nucle- 
ons. In comparison with the exclusive (e, e'p) process, the 
inclusive (e, e') scattering corresponds to an integral over 
all available nuclear states and consequently provides less 
specific information, but it is more directly related to 
the dynamics of the initial nuclear ground state. The 
width of the QE peak can give a direct measurement of 
the average momentum of nucleons in nuclei, the shape 
depends on the distribution in energy and momentum 
of the initially bound nucleons. Precise measurements 
can give direct access to integrated properties of the nu- 
clear spectral function which describes this distribution. 
A considerably body of QE data for light-to-heavy nu- 
clei in different kinematic situations has been collected 
[5, 6, 18]. Not only differential cross sections, but also 
the contribution of the separate longitudinal and trans- 
verse response functions have been considered. From the 
theoretical point of view, many efforts have been devoted 
to the description of the available data and important 
progress has been achieved in terms of experimental re- 
sults and theoretical understanding [4-6]. 

The use of the electron probe can be extended to exotic 
nuclei. The detailed study of the properties of nuclei 
far from the stability line and the evolution of nuclear 
properties with respect to the asymmetry between the 
number of neutrons and protons is one of the major topics 
of interest in modern nuclear physics. 

In the next years the advent of radioactive ion beams 
(RIB) facilities [19-21] will provide a large amount of 
data on unstable nuclei. A new generation of electron- 
RIB colliders using storage rings is under construction at 
RIKEN (Japan) [22-24] and GSI (Germany) [25]. These 
facilities will offer unprecedented opportunities to study 
the structure of exotic unstable nuclei through electron 
scattering in the ELISc experiment at FAIR in Ger- 
many [26-28] and the SCRIT project in Japan [29, 30]. 
Kincmatically complete experiments, where, in contrast 
to conventional electron scattering, all target-like reac- 
tion products are detected, will become feasible for the 
first time, allowing a clean separation of different reac- 
tion channels as well as a reduction of the unavoidable 
radiative background seen in conventional experiments. 
Therefore, even applications using stable isotope beams 
will be of interest. 



Several papers devoted to theoretical treatments of 
electron scattering off exotic nuclei have recently been 
published, underlining the usefulness of electron scatter- 
ing for investigating the structure of unstable nuclei [31- 
48]. In this work we give another contribution to this 
field. Here we present and discuss numerical predictions 
for clastic and inclusive QE electron scattering cross sec- 
tions on oxygen (^^"^^0) and calcium ("^^"^^Ca) isotopic 
chains. 

The study of the evolution of nuclear properties along 
isotopic chains requires a good knowledge of nuclear mat- 
ter distributions for protons and neutrons separately. 
The ground state densities reflect the basic properties 
of effective nuclear forces and provide fundamental nu- 
clear structure information. Elastic electron scattering 
allows to measure with excellent precision only charge 
densities and therefore proton distributions. It is much 
more difficult to measure neutron distributions. Our 
present knowledge of neutron densities comes primar- 
ily from hadron scattering experiments, the analysis 
of which requires always model-dependent assumptions 
about strong nuclear forces at low-energies. A model- 
independent probe of neutron densities is provided by 
parity-violating elastic electron scattering, where direct 
information on the neutron density can be obtained from 
the measurement of the parity- violating asymmetry Apy 
parameter, which is defined as the difference between the 
cross sections for the scattering of right- and left-handed 
longitudinally polarized electrons [49-51]. This quantity 
is related to the radius of the neutron distribution R„, 
because Z'^-boson exchange, which mediates the weak 
neutral interaction, couples mainly to neutrons and pro- 
vides a robust model- independent measurement of i?„. 

In 2012, the first measurement of Ap^ [52] (and weak 
charge form factor [53]) in the elastic scattering of po- 
larized electrons from ^''^Pb has been performed in Hall 
A at the Jefferson Lab (experiment PREX). The PREX 
collaboration obtained Apy = 0.656 ± 0.060(stat) ± 
0.014(syst) ppm, corresponding to a difference between 
the radii of the neutron and proton distributions i?„ — Rp 
= 0.33 jio^is fm. Unfortunately with such large uncer- 
tainties it is not possible to draw definite conclusions 
about the radius and the distribution of neutrons in a 
heavy finite nucleus like ^°®Pb. The problems that af- 
fected the original setup will be strongly reduced by im- 
proving electronics and radiation protection (see the re- 
cently approved PREX-II experiment [54]). In addition, 
the CREX experiment [55], with the goal of measuring 
the neutron skin of ''^Ca, has also been conditionally ap- 
proved. The combined analysis of the PREX and CREX 
experiments will allow to test the assumptions of micro- 
scopic models by testing the dependence of i?„ on the 
atomic mass number A. 

From a theoretical point of view, parity violation in 
elastic electron scattering has been recently studied by 
several authors within the framework of mean field ap- 
proaches [43-46, 56-59] after almost a decade from the 
first detailed calculations [51, 60, 61]. In particular, in 
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Ref. [58] the authors suggested that a 1% measurement 
of Ap^ can constrain the slope L of the symmetry energy 
close to a 10 MeV level of accuracy. A precise measure- 
ment of Rn is definitely a crucial step to improve our 
knowledge of neutron-rich matter, i.e., the outer part of 
neutron stars. 

In this paper we present calculations of Apy performed 
for oxygen and calcium isotopic chains and test the iso- 
topic dependence. In addition, we present a comparison 
with the results of the PREX collaboration and provide 
some estimates for the future experiment CREX. 

The basic ingredients of the calculations for both elas- 
tic and QE scattering are the ground state wave functions 
of proton and neutron s.p. states, for which we have used 
a relativistic mean field (RMF) model. In the last years 
mean field approaches have been very successfully em- 
ployed in several aspects of nuclear structure phenom- 
ena, in particular when density-dependent parameteri- 
zations are considered. Effective hadron field theories 
with medium dependent meson-nuclcon vertices retain 
the basic structure of the relativistic mean-field frame- 
work, but can be more directly related to the underlying 
microscopic description of nuclear interactions [62-64]. 
After the first applications [65, 66], restricted to infinite 
systems and spherical nuclei, with the parametcrizations 
DD-MEl [67] and DD-ME2 [68], calculations have been 
extended to open shell nuclei, exotic systems, superheav- 
ies, and collective resonances [69]. 

The cross sections for elastic electron scattering are 
obtained by summing up the phase shifts which result 
from the numerical solution of the partial wave Dirac 
equation and includes Coulomb distortion effects. 

In inclusive QE electron scattering a proper description 
of the final-state interactions (FSI) between the emitted 
nucleon and the residual nucleus is an essential ingredi- 
ent for the comparison with data. For the calculations 
presented in this paper we have employed the relativis- 
tic Green's function model, which has been widely and 
successfully applied to the analysis of QE electron and 
neutrino-nucleus scattering data on different nuclei [70- 
77]. 

As a first step, we test the results of the models in com- 
parison with some of the data available for elastic and QE 
scattering on ^^O and '*°Ca nuclei. Then, with extensive 
calculations on oxygen and calcium isotopic chains, we 
investigate the evolution of some ground-state properties 
as a function of the neutron number. 

The case of the exclusive quasi-free (e, e'p) reaction has 
been investigated in [47, 78] , where the cross sections ob- 
tained with different relativistic and nonrelativistic ap- 
proaches based on the mean-field description for the pro- 
ton bound state wave function are compared for oxygen 
and calcium isotopic chains. 

The paper is organized as follows. In Sec. II we give 
the basic formalism involved in the description of elastic, 
quasi-elastic, and parity-violating electron scattering. In 
Sec. Ill we present a brief discussion of the relativis- 
tic mean-field model and of the calculations of the self- 



consistent ground-state proton and neutron densities of 
calcium and oxygen isotopes. In Sec. IV we outline the 
main features of the relativistic Green's function model, 
which is used to describe final-state interactions in the 
inclusive quasi-elastic electron scattering. In Sec. V we 
show and discuss our theoretical results obtained for elas- 
tic, quasi-clastic, and parity- violating electron scattering. 
Finally, in Sec. VI we summarize our results and present 
our conclusions. 



II. ELASTIC AND QUASI-ELASTIC 
ELECTRON SCATTERING 

A. Elastic electron scattering 

In the one-photon exchange approximation and ne- 
glecting the effect of the nuclear Coulomb field on in- 
coming and outgoing electrons, i.e. in the plane- wave 
Born approximation (PWBA), the differential cross sec- 
tion for the clastic scattering of an electron off a spherical 
spin-zero nucleus is given by 

where ft' is the scattered electron solid angle, am is the 
Mott cross section [4, 5] and 

FpiQ') = ^ J dr j,{qr)pp{r) , (2) 

is the charge form factor for a spherical nuclear charge 
(point proton) density Pp{r), = |q|^ — w^, and = 
(a;, q) is the four-momentum transfer. 

The PWBA is, however, not adequate for medium and 
heavy nuclei where the distortion produced on the elec- 
tron wave functions by the nuclear Coulomb potential 
V{r) from Pp{r) can have significant effects. The DWBA 
cross sections are obtained by summing up the phase 
shifts which result from the numerical solutions of the 
partial wave Dirac equation. 

B. Inclusive quasi-elastic electron scattering 

In the one-photon exchange approximation the inclu- 
sive differential cross section for the QE (e, e') scattering 
on a nucleus is obtained from the contraction between 
the lepton and hadron tensors as [5] 

{-J^)^^--M[v,R, + vrRr] , (3) 

where e' is the energy of the scattered electron. The coef- 
ficients V come from the components of the lepton tensor 
that, under the assumption of the plane-wave approxi- 
mation for the electron wave functions, depend only on 
the lepton kinematics. 
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where 9 is the electron scattering angle. All nuclear 
structure information is contained in the longitudinal and 
transverse response functions Rl and Rt, expressed by 



RL{q,u:)^W'''\q,uj) , 

RT{q.u:) = W^\q,uj) + W^''{q,uj) , 



(5) 



in terms of the diagonal components of the hadron tensor, 
that is given by bilinear products of the transition matrix 
elements of the nuclear electromagnetic many-body cur- 
rent operator between the initial state of the nucleus 
I V^o), of energy Eq, and the final states | ^/), of energy 
Ef, both eigenstates of the nuclear Hamiltonian H, as 



i f 

X (vfo I J^Hq) I */> KEo+^-Ef), (6) 

involving an average over the initial states and a sum over 
the undetected final states. The sum runs over the scat- 
tering states corresponding to all of the allowed asymp- 
totic configurations and includes possible discrete states 
[4, 5]. 

The hadron tensor can cquivalently be expressed as 

= -ilm(vI/o I J''Hq)G{Ef)J^^{q) \ ^o) (7) 

where Ef = Eq+u and G{Ef) is the many-body Green's 
function related to the many-body nuclear Hamiltonian 
H. 

The hadron in tensor in Eq. (7) contains the full many- 
body propagator of the nuclear system. As such, it is an 
extremely complicated object and some approximations 
are needed to reduce the calculation of the nuclear re- 
sponse to a tractable form. 



C. Parity-violating electron scattering 

When a photon is exchanged between two charged par- 
ticles a Z° boson is also exchanged. At the energies of 
interest in electron scattering the strength of the weak 
process mediated by the Z° boson is negligible compared 
with the electromagnetic strength. The role played by 
the Z'^ exchange is therefore not significant unless an ex- 
periment is set up to measure a parity-violating observ- 
able. While the electromagnetic interaction conserves 
parity, the weak interaction does not and this is how we 
are sensitive to Z° exchange in electron scattering. 

The degree of parity violation can be measured by the 
parity- violating asymmetry Apy^ or helicity asymmetry, 
which is defined as the difference between the cross sec- 
tions for the scattering of electrons longitudinally polar- 
ized parallel and antiparallel to their momentum. This 
difference arises from the interference of photon and Z'^ 
exchange. As it has been shown in Refs. [49, 60], the 
asymmetry in the parity violating elastic polarized elec- 
tron scattering represents an almost direct measurement 



of the Fourier transform of the neutron density. Follow- 
ing Ref. [51], we consider the elastic electron scattering 
on a spin-zero nucleus, i.e. on the potential 



U{r) = V{r)+-f5A{r), 



(8) 



where V{r) is the Coulomb potential, and A{r) results 
from the weak neutral current amplitude 



A{r) = Gi./23/Vw(r), 



(9) 



where Gp ^ 1.16639 x 10"^^ MeV"^ jg ^he Fermi con- 
stant. 

The weak charge density is defined 

pw{r)=j dr' Ge {\r-r'\) 

X [-Pn{r') + (1 - Asm^ew)Pp{r')] (10) 

where pn and pp are point neutron and proton densi- 
ties. The electric form factor of the proton is Gsir) ~ 
A3/(87r)e-^^ with A = 4.27 fni"!, and sin^e^ ^ 0.23 
for the Weinberg angle. The weak potential is therefore 
dominated by the neutron density and this is why parity- 
violating electron scattering (PVES) has been suggested 
to study neutron densities in nuclei. 

In the limit of vanishing electron mass, the helicity 
states 'ii± = (1/2)(1 ± 75)^* satisfy the Dirac equation 



with 



[a .p + C/±(r)]*± =-B^'±, 



U±{r) = V{r)±A{r). 



(11) 



(12) 



The parity- violating asymmetry Apy, or helicity asym- 
metry, is defined 



da+ 1 d£t — da- / dil 
^ da+/dn + da-/dn 



(13) 



where -l-(— ) refers to the clastic scattering on the poten- 
tial U±{r). 

In QE electron scattering the helicity asymmetry is 
obtained in terms of kinematic coefficients v from the 
lepton tensor and of nuclear response functions as [72] 



„, TjAV I „, dAV I idVA 

^pv = ^0 5 — ; 5 



The factor Aq is defined as 



^0 



Gf 
2\/2- 



7ra 



(14) 



(15) 



where a is the fine structure constant. The denominator 
in Eq. (14) contains the parity-conserving cross section of 
Eq. (3), the numerator the parity- violating contribution. 



tan ■ 



'tan^ 



(16) 
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and the response functions R are given in terms of the 
polarized components of the hadron tensor W^'^ [72] as 



tdVA 

lirp 



11 



igv - WD 



(17) 



where the superscript AV denotes interference of axial- 
vector leptonic current with vector hadronic current (the 
reverse for VA) and the couplings gA = —1/2 and gv = 
-l/2 + 2sin2eiv ~ -0.04. 



III. RELATIVISTIC MODEL FOR 
GROUND-STATE OBSERVABLES 

In the standard representation of rclativistic mean field 
approaches the nucleus is described as a system of Dirac 
nucleons coupled to the exchange mesons and the elec- 
tromagnetic field through an effective Lagrangian. The 
isoscalar scalar-meson (a), the isoscalar vector- meson 
(w), and the isovector vector- meson (p) build the minimal 
set of meson fields that together with the electromagnetic 
field (7) is necessary for a quantitative description of bulk 
and s.p. nuclear properties. The model is defined by the 
Lagrangian density 



C = C 



N 



^7) 



(18) 



where /Zat denotes the Lagrangian of the free nucleon, 
Lra is the Lagrangian of the free meson fields and the 
simplest set of interaction terms is contained in Hint '■ 

Cint = -gai'o-ip - gu^ipj'^ujf^ip - gpipj'^f ■ Pf_,ip . (19) 

The couplings of the cr-meson and w-meson to the nucleon 
are assumed to be of the form: 



9i{p) = 9i{Psat)fi{x) for i 



where 



l + h{x + dif 



(20) 



(21) 



is a function oi x = p/ psat, and psat denotes the nucleon 
density at saturation in symmetric nuclear matter. Con- 
straints at nuclear matter saturation density and at zero 
density are used to reduce the number of independent 
parameters in Eq. (21) to three. Three additional pa- 
rameters in the isoscalar channel are gaipsat), gcjiPsat), 
and nia , that is the mass of the phenomenological a me- 
son. For the p meson coupling the functional form of 
the density dependence is suggested by Dirac-Brueckner 
calculations of asymmetric nuclear matter: 



9p{p) 9p{Psat) exp[-ap(a; - 1)] 



(22) 



and the isovector channel is parametrized by gp{psat) and 
a p. Bare values are used for the masses of the to and p 



mesons: m,. 



is determined by eight independent parameters, adjusted 
to the properties of symmetric and asymmetric nuclear 
matter, binding energies, charge radii, and neutron radii 
of spherical nuclei [68]. The interaction has been tested 
in the calculation of ground state properties of a large set 
of spherical and deformed nuclei. When used in the rcl- 
ativistic random-phase approximation, DD-ME2 repro- 
duces with high accuracy data on isoscalar and isovector 
collective excitations. For open-shell nuclei we employed 
a schematic ansatz: the constant gap approximation with 
empirical A given by the 5-point formula [79] 



A(5)(A^o) 



E{Na 



AE{No + 1) + &E{Nq) 



AE{No - 1) + E{Na - 2) 



(23) 



783 MeV and = 763 MeV. DD-ME2 



In Figs. 1 and 2 we plot the neutron (proton) density dis- 
tributions Pn(p) as a function of the radial coordinate r 
for oxygen and calcium isotopes, respectively. These den- 
sity distributions are the sum of the squared moduli of 
the s.p. neutron (proton) wave functions. All the nuclei 
we have investigated resulted to be bound. From the ex- 
perimental point of view it seems rather well established 
that the neutron drip line for the oxygen isotopes starts 
with ^^O [80] and, therefore, ^^O should not be bound. 

For both isotopic chains the differences between the 
proton and neutron densities in Figs. 1 and 2 generally 
increase with the neutron number. When the number 
of neutrons increases there is a gradual increase of the 
neutron radius. The differences of the neutron density 
profiles in the nuclear interior display pronounced shell 
effects. The effect of adding neutrons is to populate and 
extend the neutron densities and, to a minor extent, also 
the proton densities. In the case of protons, however, 
there is a decrease of the density in the nuclear interior 
to preserve the normalization to the constant number of 
protons. 



IV. RELATIVISTIC GREEN'S FUNCTION 
MODEL FOR INCLUSIVE QUASI-ELASTIC 
ELECTRON SCATTERING 

In the QE region the nuclear response is dominated by 
one-nucleon knockout processes, where the scattering oc- 
curs with only one nucleon that is subsequently emitted. 
The remaining nucleons of the target behave as simple 
spectators and QE electron scattering can adequately be 
described in the rclativistic impulse approximation (RIA) 
by the sum of incoherent processes involving only one nu- 
cleon scattering and the components of the hadron ten- 
sor of Eq. (6) are obtained from the sum, over all the s.p. 
shell-model states, of the squared absolute value of the 
transition matrix elements of the single-nucleon current. 

A reliable description of final-state interactions (FSI) 
between the ejected nucleon and the residual nucleus is 
an essential ingredient for the comparison with data. 
In the case of exclusive (e, e'p) processes, the use of 
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complex optical potentials in the rclativistic distorted 
wave impulse approximation (RDWIA) has been able 
to successfully describe a wide number of experimental 
data [16, 47, 48, 81-85]. It is clear that the pure RDWIA 
approach, based on the use of an absorptive complex po- 
tential, would be inconsistent in the analysis of inclusive 
scattering, where all final-state channels should be re- 
tained and the total flux, although redistributed among 
all possible channels due to FSI, must be conserved. Dif- 
ferent approaches have been used to describe FSI in RIA 
calculations for the inclusive QE electron- and neutrino- 
nucleus scattering [70-75, 77, 86-95]. In the rclativis- 
tic plane- wave impulse approximation (RP WIA) , FSI are 
simply neglected. In another approach, FSI are included 
in calculations where the final nucleon state is evaluated 
with real potentials, cither retaining only the real part of 
the rclativistic energy-dependent complex optical poten- 
tial (rROP) [94, 95], or using the same rclativistic mean 
field potential considered in describing the initial nucleon 
state [90]. 

In the rclativistic Green's function (RGF) model FSI 
arc described in the inclusive process consistently with 
the exclusive scattering by the same complex optical po- 
tential, but the imaginary part is used in the two cases 
in a different way and in the inclusive scattering it redis- 
tributes the flux in all the channels and the total flux is 
conserved. Detailed discussions of the RGF model can 
be found in Refs. [70-77, 96, 97]. 

In the RGF model with suitable approximations, which 
are mainly related to the impulse approximation, the 
components of the nuclear response of Eq. (7) are writ- 
ten in terms of the s.p. optical model Green's function. 
The spectral representation of the s.p. Green's function, 
which is based on a biorthogonal expansion in terms of a 
non-Hermitian optical potential and of its Hermitian con- 
jugate, can be exploited to avoid the explicit calculation 
of the s.p. Green's function and obtain the components 
of the hadron tensor in the form [70] 



Re Til^^iEf 



1 

-V 

IT 



M 



E 



f 



E 



Im TZi£,Ef-en) 



(?4) 



where V denotes the principal value of the integral, n is 
the cigenstatc of the residual nucleus with energy £„, and 



TnS.E) ^ A„(^„ I j^Hq)^l~ViE) I x^iE)) 

x{xi-\E) I ^l-V{E)f{q) I ^,,) (25) 

The factor ^1 - V'{E), where V'iE) is the energy 
derivative of the optical potential, accounts for interfer- 
ence effects between different channels and justifies the 
replacement in the calculations of the Feshbach optical 
potential V of the RGF model, for which neither mi- 
croscopic nor empirical calculations are available, by the 
local phenomenological optical potential [70, 96, 97]. 



Disregarding the square root correction, the second 
matrix element in Eq. (25) is the transition amplitude 
of the exclusive single-nucleon knockout, which appears 
in usual RDWIA calculations. In this matrix element 
is the one-body nuclear current, x*-"-* is the s.p. scat- 
tering state of the emitted nucleon with energy £, tpn is 
the overlap between the ground state of the target and 
the final state ri, i.e., a s.p. bound state, and the spec- 
troscopic factor A„ is the norm of the overlap function. 
The s.p. bound and scattering states, ipn and x^~\ 
eigenfunctions of a non Hermitian energy- dependent op- 
tical potential and of its Hermitian conjugate at different 
energies. 

In the exclusive one-nuclcon knockout the imaginary 
part of the optical potential accounts for the flux lost 
in the channel n towards the channels different from n, 
which are not included in the exclusive process. In the in- 
clusive response, where all the channels are included, this 
loss is compensated by a corresponding gain of flux due 
to the flux lost, towards the channel n, in the other final 
states asymptotically originated by the channels differ- 
ent from n. This compensation is performed by the first 
matrix element in the right hand side of Eq. (25), which 
involves the eigenfunction Xs (E) of the Hermitian con- 
jugate optical potential, where the imaginary part has an 
opposite sign and has the effect of increasing the strength. 
Therefore, in the RGF approach the imaginary part of 
the optical potential redistributes the flux lost in a chan- 
nel in the other channels, and in the sum over n the total 
flux is conserved. If the imaginary part of the optical po- 
tential is neglected, the second term in Eq. (24) vanishes 
and, but for the square root factor, the first term gives 
the rROP approach. 

The RGF model has been applied to parity-violating 
QE electron scattering in Ref. [72]. The main steps 
of the model are the same, the expressions for the 
electromagnetic-weak interference components of the 
hadron tensor Wj'^ can be found in Ref. [72] . 



V. RESULTS 

In this section we present and discuss numerical pre- 
dictions for elastic and QE electron scattering which 
can hopefully be useful for future measurements in ex- 
perimental RIB facilities. We study the evolution of 
some electron scattering observables in isotopic chains 
of medium systems, which are exemplified by the cases 
of the oxygen and calcium isotopes. Many of these nuclei 
lie in the region of the nuclear chart that is likely to be 
explored in future electron-scattering experiments. As 
a first step, with a few numerical examples we test the 
results of our models in comparison with available data. 
Then, for each nucleus in an isotopic chain, we compute 
and compare the associated elastic and QE cross sections 
and parity- violating asymmetry in order to obtain infor- 
mation on the effects of isospin asymmetry on nuclear 
structure. 
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A. Elastic electron scattering 

The cross sections for clastic electron scattering have 
been calculated in the DWBA and with the self- 
consistent relativistic ground state charge densities de- 
scribed in Sect. III. In the PWBA the cross section 
is proportional to the Fourier transform of the proton 
charge density (see Eq. (2)) and reflects its behavior also 
when Coulomb distortion is included in the calculations. 
In different studies of the charge form factors along iso- 
topic chains [36, 40, 41, 98, 99] it has been found that, 
when the number of neutrons increases, the squared mod- 
ulus of the charge form factor and the position of its min- 
ima show, respectively, an upward trend and a significant 
inward shifting in the momentum transfer. 

An example of the comparison between theoretical and 
experimental differential cross sections is displayed in 
Fig. 3 for clastic electron scattering on -'^^O at an electron 
energy e = 374.5 MeV and on ""Ca at e = 496.8 MeV. 
The general trend of the experimental data is reasonably 
reproduced by the calculations. Both experimental cross 
sections considered in the figure are well described at low 
scattering angles. For ""^Ca there is a fair agreement be- 
tween theory and data also at larger angles, while for ^^O 
data beyond the minimum are somewhat underestimated 
by the theoretical results. 

The calculated differential cross sections for elastic 
electron scattering on various oxygen isotopes (^^^^^O) 
at £ = 374.5 MeV and on calcium isotopes ('^^"^^Ca) at 
e = 496.8.5 MeV are shown in Figs. 4 and 5, respec- 
tively. With increasing neutron number the positions 
of the diffraction minima shift toward smaller scatter- 
ing angles, i.e., towards smaller values of the momentum 
transfer. The shift indicates that increasing the neutron 
number the skin charge densities increase and the charge 
density distributions tend to be less extended [40]. The 
shift of the minima towards smaller q is in general ac- 
companied by a simultaneous increase in the height of 
the maxima. The behavior is similar for both isotopic 
chains here considered and is in agreement with the re- 
sults found in previous studies of charge form factors on 
various isotopic chains, which were carried out with dif- 
ferent mean-field models [36, 40, 41, 98, 99]. 

B. Quasi-elastic electron scattering 

The cross sections for QE electron scattering have been 
computed with the RGF model. Some results obtained in 
the RPWIA arc also presented for a comparison. In the 
calculations of the matrix elements in Eq. (25) the s.p. 
bound nuclcon states are obtained from the relativistic 
mean-field model with density-dependent meson-nucleon 
vertices and the DD-ME2 parametrization as described 
in Sect. III. The s.p. scattering states are eigenfunc- 
tions of the energy-dependent and A— dependent [A is 
the mass number) parameterization for the relativistic 
optical potential of Rcf. [100], which is fitted to pro- 



ton elastic scattering data on several nuclei in an energy 
range up to 1040 MeV. For the single- nucleon current we 
have used the relativistic free nucleon expression denoted 
as CC2 [70, 101]. 

The predictions of the RGF model have been compared 
with experimental data for QE electron- and neutrino- 
nucleus scattering in a series of papers [70-77, 88, 92], 
where the calculations have been performed with differ- 
ent relativistic mean-field models for the bound states 
and different paramcterizations of the relativistic optical 
potential. 

In Fig. 6 our RGF results are compared with the exper- 
imental (e, e') cross sections for two different kinematics 
on ^''O and ^°Ca target nuclei [102, 103]. The agreement 
with the data is satisfactory, at least in the energy re- 
gion of the QE peak. The RGF model was developed 
to describe FSI in inclusive QE electron scattering and 
is in general able to give a reasonable and even good 
description of QE data. For energy regions below and 
above the QE peak other contributions, not included in 
the RGF model, can be important. Even in the QE re- 
gion, the relevance of contributions like meson exchange 
currents and Delta effects should be carefully evaluated 
before definite conclusions can be drawn about the com- 
parison with data [104-106]. Such contributions may be 
significant even in the QE region, in particular in kine- 
matics where the transverse component of the nuclear 
response plays a major role in the cross section. 

The cross section of the inclusive QE (e, e') reaction on 
"-280 isotopes at e 1080 MeV and 9 = 32° are shown 
in Fig. 7. In a first approximation, we have neglected FSI 
and calculations have been performed in the RPWIA. In 
this case, the differences between the results for the var- 
ious isotopes are entirely due to the differences in the 
s.p. bound state wave functions of each isotope. While 
only the charge proton density distribution contributes to 
the cross section of elastic electron scattering, the cross 
section of QE electron scattering is obtained from the 
sum of all the integrated exclusive one-nucleon knockout 
processes, due to the interaction of the probe with all the 
individual nucleons, protons and neutrons, of the nucleus 
and contains information on the dynamics of the initial 
nuclear ground state. The separate contributions from 
protons and neutrons are also shown in the lower panel 
of Fig. 7. In an usual experiment where only the scat- 
tered electron is detected these two quantities cannot be 
separated experimentally, but their comparison can give 
useful information on the different role of protons and 
neutrons in the inclusive QE cross section. The main 
role is played by protons, which give most of the contri- 
bution. Increasing the neutron number it is quite natural 
to understand the proportional increase of the neutron 
contribution. No significant increase is found in the pro- 
ton contribution. Thus, the increase of the cross section 
in the upper panel of the figure is due to the increase of 
the neutron contribution. The shift of the proton con- 
tribution towards higher values of uj seen in the figure 
is mainly related to the increase of the proton separa- 
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tion energy with increasing neutron number (increasing 
the neutron number the protons experience more binding 
and their separation energies increase) than to changes 
in the proton wave functions. A different and opposite 
shift can be seen in the case of the neutron contribution 
and, therefore, the final effect is that the shift is strongly 
reduced in the QE cross section shown in the upper panel 
of Fig. 7. 

In Fig. 8 we show the QE (e, e') cross sections calcu- 
lated for oxygen isotopes with the RGF model and in the 
same kinematics as in Fig. 7. The general trend of the 
cross sections, their magnitude, and their evolution with 
respect to the change of the neutron number are gen- 
erally similar in RPWIA and RGF. The FSI effects in 
the RGF calculations produce, however, some differences 
which can be seen in the low energy transferred region, 
where the cross sections for i4,i6,18q ^^.^ enhanced with 
respect to those for ^^O and ^^O. In addition, the shift 
towards higher uj is more significant than in the RPWIA 
case, but for ^®0. 

The cross section of the inclusive QE (e, e') reaction 
on 36-56^3^ isotopes at e = 560 MeV and 6 = 60° cal- 
culated in the RPWIA and in the RGF are shown in 
Figs. 9 and 10, respectively. The general behavior of the 
cross sections and their evolution with increasing neutron 
number is similar for calcium and oxygen isotopes. The 
magnitude increases with the neutron number, but FSI 
effects are somewhat more visible for calcium isotopes. 
The RGF cross sections on 36,40,44,48Qg^ 
enlarged over a wide range of w and are slightly reduced 
with respect to the RPWIA results in Fig. 9. This is 
particularly visible for ^^Ca and produces an apparently 
large gap between the cross sections of '^^Ca and ^^Ca. 

As a final comment, we can add that interesting and 
peculiar effects are obtained in the evolution of QE in- 
clusive cross section along isotopic chains, but it is not 
easy to relate them to changes in the matter distribution, 
which can be significant, particularly in the center of the 
nucleus. 



C. Parity-violating asymmetry 

The calculation starts with the self-consistent relativis- 
tic ground state proton and neutron densities (see Sect. 
III). The charge and weak densities are calculated by 
folding the point proton and neutron densities (see Eq. 
(10)). The resulting Coulomb potential V{r) and weak 
potential A{r) (sec Eq.(9)) are used to construct U±{r). 
The cross sections for elastic electron scattering are ob- 
tained by summing up the phase shifts which result from 
the numerical solution of the partial wave Dirac equation. 
The calculation includes Coulomb distortion effects. The 
cross sections for positive and negative helicity electron 
states are calculated, and the resulting asymmetry pa- 
rameter Apy is plotted as a function of the scattering 
angle. 

In Figs. 11 and 12 wc plot the parity-violating asym- 



metry parameters Apy for 14-28 q ^^^^ 36-56(-;j^ nuclei for 
elastic electron scattering at e — 850 MeV. At e = 850 
MeV the values of Apy are of the order of 10^^ for rela- 
tively small angles. 

As suggested in Ref. [49] the asymmetry parameter 
Apy provides a direct measurement of the Fourier trans- 
form of the neutron density. This relation has been tested 
and confirmed in Ref. [60] comparing asymmetries and 
the squares of the Fourier transforms of the neutron den- 
sities. Another way to relate Apy to neutron distributions 
of finite nuclei is by looking at possible linear correlations 
between the asymmetry parameter and some well defined 
observables. We suggest to use the first minima positions 
9,nin and the neutron excess A = {N—Z)/Z, i.e., how the 
minima of Apy evolve from neutron-poor to neutron-rich 
nuclei. In fact, by neglecting strangeness contributions 
and the electric neutron form factor, the parity violating 
asymmetry can be rewritten as 



An 



4sin^ ew + ^ + 

Zj 



(26) 



where k = Gf / ['^i^oiV^), and R contains nuclear struc- 
ture effects through neutron/proton form factors R = 
F°{q)/F^{q) - 1 [107]. In Fig. 13, 0™„ is plotted as 
a function of A for oxygen and calcium isotopes. The 
dashed lines suggest that for both isotope chains the evo- 
lution of Apy as function of A is well approximated by a 
linear fit with a very similar slope. To test the robust- 
ness of this correlation it is interesting to study if Apy 
is affected by density distribution oscillations at small 
radii that could appear in some selected cases. In Ref. 
[39], ^^O and ^'*0 isotopes have been studied as possible 
candidates for "bubble" nuclei, i.e., nuclear systems with 
a strong depicted central density. In Fig. 14 we plot 
the asymmetry parameters Apy for these nuclei. Neu- 
tron density profiles show large differences at small dis- 
tances but with no appreciable effects in the correspond- 
ing asymmetries. As a consequence, we can safely assume 
that Apy is a reliable observable to study neutron radii. 
The inclusion of pairing correlations does not modify this 
assertion. 

In addition to predictions about oxygen and calcium 
isotope chains wc also provide calculations for recent 
measurements and future experiments. In Fig. 15 we 
show our theoretical predictions for the empirical values 
extracted from the first run of the PREX experiment on 
208pb at e = 1.06 GeV. In Ref. [53] the weak charge 
density {—pw) has been deduced from the weak charge 
form factor. The error band (shaded area) represents the 
incoherent sum of experimental and model errors. Our 
prediction, plotted by the red line in the left panel, is 
in rather good agreement with empirical data. In fact, 
if we evaluate the corresponding asymmetry parameter 
Apy averaged over the acceptance function e{0) [108] 



{Apy) — 



jA6 sine Apy{e)%e{e) 



SAO sin 0^6(0) 



(27) 
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we find 0.712 ppm, in very good agreement with the em- 
pirical estimate 0.656 ± 0.060(stat) ± 0.014(syst) ppm. 
In Fig. 16 we calculate the asymmetry parameter for 
''^Ca with 2.2 GeV electrons as planned for the CREX 
experiment. For energies well above the 1 GcV region, of 
course, the elastic scattering approximation is not com- 
pletely under control and corrections due to possible in- 
elasticities should be taken into account. We plan to 
extend our calculation in a forthcoming paper. 

As an example of the parity-violating asymmetry 
for QE scattering, in Fig. 17 we show Apy for the 
14,16,18,22,280 isotopes evaluated with the RGF in the 
same kinematics as in Fig. 8. Note that the results are 
rescaled by the factor 10^. The RPWIA results are al- 
ways similar to the RGF ones and are not presented here. 
The asymmetry is almost constant to a few x 10~^ for 
^"^O and ^^O; whereas for ^^O and ^^0 it goes up to 
w —2 X 10^^ in the low energy transferred region. There 
are visible relative differences between the asymmetries 
for the oxygen isotopes that we have considered, but we 
are aware that a measurement of Ap^, in QE electron scat- 
tering off finite nuclei is extremely challenging with the 
presently available facilities. 



VI. SUMMARY AND CONCLUSIONS 

We have presented and discussed numerical predictions 
for the cross section and the parity-violating asymmetry 
in elastic and quasi-elastic electron scattering on oxygen 
and calcium isotope chains with the aim to investigate 
their evolution with increasing neutron number. 

The understanding of the properties of exotic nuclei 
is one of the major topic of interest in modern nuclear 
physics. Large efforts in this directions have been done 
over last years and are planned for the future. The use 
of electrons as probe provides a powerful tool to achieve 
this goal. The RIB facilities in different laboratories have 
opened the possibility to give insight into nuclear struc- 
tures which are not available in nature, as they are not 
stable, but which are important in astrophysics and had 
a relevant role in the nucleosynthesis. 

Electron scattering is well fitted for studying nuclear 
properties, as its interaction is well known and relatively 
weak with respect to the hadron force and can there- 
fore more adequately explore the details of inner nu- 
clear structures. As a consequence of this weakness, the 
cross sections become very small and more difficult ex- 
periments have to be performed. Electron scattering ex- 
periments off exotic nuclei have been proposed in the 
ELISe experiment at FAIR and in the SCRIT project at 
RIKEN. We hope that the existing proposals will be con- 
sidered and approved in next years. Our theoretical pre- 
dictions will be useful for clarifying the different aspects 
of the measurements, giving information on the order 
of magnitude of the measurable quantities and therefore 
making possible a more precise evaluation of the experi- 
mental difficulties. Moreover, a theoretical investigation 



can be helpful to envisage the most interesting quanti- 
ties to be measured in order to explore the properties of 
exotic nuclear structures. 

In this work, both elastic and inclusive quasi-elastic 
electron scattering have been considered. The elastic 
scattering can give information on the global properties 
of nuclei and, in particular, on the different behavior of 
proton and neutron density distributions. The inclusive 
quasi-elastic scattering is affected by the dynamical prop- 
erties, being the integral of the spectral density function 
over all the available final states, and, due to the reaction 
mechanism, preferably exploits the single particle aspects 
of the nucleus. In addition, when combined with the ex- 
clusive (e, e'p) scattering, it is able to explore the evolu- 
tion of the single particle model with increasing asymme- 
try between the number of neutrons and protons. Many 
interesting phenomena are predicted in this situation: 
in particular, the modification of the shell model magic 
numbers. A definite response can be obtained from the 
comparison with experimental data, which will discrim- 
inate between the different theoretical models, mainly 
referring to rclativistic mean field approaches. 

As case studies for the present investigation we have 
selected oxygen and calcium isotope chains. The cal- 
culations have been carried out within the framework 
of the relativistic mean field model. The nuclear wave 
functions are obtained considering a system of nucleons 
coupled to the exchange mesons and the electromagnetic 
field through an effective Lagrangian. The calculated 
cross sections include both the hadronic and Coulomb 
final states interactions. The inclusive quasi-elastic scat- 
tering is calculated with the relativistic Green's function 
model, which conserves the global particle fiux in all the 
final state channels, as it is required in an inclusive reac- 
tion. 

First, the models have been compared with experimen- 
tal data already available on stable isotopes in order to 
check their reliability. Then, the same models have been 
used to calculate elastic and inclusive quasi-elastic cross 
sections on exotic isotopes chains. The possible disagree- 
ment of the experimental findings from the theoretical 
predictions will be a clear indication of the insurgence of 
new phenomena related to the proton to neutron asym- 
metry. 

Our results show an evolution of the calculated quanti- 
ties without discontinuities. The increase of the neutron 
number essentially produces an increase of the nuclear 
and proton densities and a flattening of the charge den- 
sity. 

The parity-violating asymmetry parameter has been 
calculated in order to investigate the neutron skin, as 
the weak current is essentially obtained from the inter- 
action with neutrons. Numerical predictions have been 
presented and discussed for oxygen and calcium isotopes. 
Moreover, we have compared our calculations with the 
results of the first measurement of the asymmetry pa- 
rameter achieved by the PREX Collaboration on ^°^Pb 
and have obtained a good agreement with the empirical 
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value. A prediction for the future experiment CREX on 
^*Ca has also been given. 



ACKNOWLEDGMENTS 

This work was partially supported by the Italian MIUR 
through the PRIN 2009 research project. 



[1] R. Hofstadter, Rev. Mod. Phys. 28, 214 (1956). 

[2] T. W. Donnelly and J. D. Walecka, Ann. Rev. Nucl. 

Part. Sci. 25, 329 (1975). 
[3] T. W. Donnelly and I. Sick, Rev. Mod. Phys. 56, 461 

(1984). 

[4] S. Boffi, C. Giusti, and F. D. Pacati, Phys. Rept. 226, 
1 (1993). 

[5] S. Boffi, C. Giusti, F. D. Pacati, and M. Radici, Elec- 
tromagnetic Response of Atomic Nuclei, Oxford Studies 
in Nuclear Physics, Vol. 20 (Clarendon Press, Oxford, 
1996). 

[6] O. Benhar, D. Day, and I. Sick, Rev. Mod. Phys. 80, 
189 (2008). 

[7] H. De Vries, C. W. De Jager, and C. De Vries, Atom. 

Data Nucl. Data Tabl. 36, 495 (1987). 
[8] G. Fricke, C. Bernhardt, K. Heilig, L. A. Schaller, 
L. Schellenberg, E. B. Shera, and C. W. De Jager, 
Atom. Data Nucl. Data Tabl. 60, 177 (1995). 
[9] J. Heisenberg and H. P. Blok, Ann. Rev. Nucl. Part. Sci. 
33, 569 (1983). 
[10] E. Moya de Guerra, Phys. Rept. 138, 293 (1986). 
[11] M. N. Harakeh and A. van der Woude, Giant Res- 
onances, Oxford Studies in Nuclear Physics, Vol. 24 
(Clarendon Press, Oxford, 2001). 
[12] S. FruUani and J. Mougey, Adv. Nucl. Phys. 14, 1 
(1984). 

[13] M. Bernheim, A. Bussiere, J. Mougey, D. Royer, 
S. Tarnowski, D. Turck-Chieze, S. FruUani, S. Boffi, 
C. Giusti, and F. D. Pacati, Nucl. Phys. A375, 381 
(1982). 

[14] L. Lapikas, Nucl. Phys. A553, 297 (1993). 

[15] P. K. A. de Witt Huberts, J. Phys. G16, 507 (1990). 

[16] J. M. Udi'as, P. Sarriguren, E. Moya de Guerra, E. Gar- 

rido, and J. A. Caballero, Phys. Rev. C 48, 2731 (1993), 

arXiv:nucbth/9310004 [nucl-th]. 
[17] J. Gao et al. (The Jefferson Lab Hall A Collaboration), 

Phys. Rev. Lett. 84, 3265 (2000). 
[18] http : //faculty . Virginia, edu/qes- archive/ index . 

html. 

[19] L Tanhihata, Prog. Part. Nucl. Phys. 35, 505 (1995). 

[20] H. Geissel, G. Miizenberg, and R. Riisager, Ann. Rev. 
Nucl. Part. Sci. 45, 163 (1995). 

[21] A. Mueller, Prog. Part. Nucl. Phys. 46, 359 (2001). 

[22] T. Suda, K. Maruyama, and I. Tanhihata, RIKEN Ac- 
cel. Prog. Rep. 34, 49 (2001). 

[23] T. Suda, M. Wakasugi, T. Emoto, K. Ishii, S. Ito, 
K. Kurita, A. Kuwajima, A. Noda, T. Shirai, T. Tamae, 
H. Tonga, S. Wang, and Y. Yano, Phys. Rev. Lett. 102, 
102501 (2009). 

[24] T. Katayama, T. Suda, and I. Tanhihata, Physica 
Scripta T104, 129 (2003). 



[25] "An International Accelerator Facility for Beams of Ions 
and Antiprotons" , GSI report 2006, http://www.gsi. 
de /GS I - Futur e/ cdr / . 

[26] http : //www.gsi .de/f orschung/f air_experiments/ 
elise/ index_e . html. 

[27] H. Simon, Nucl. Phys. A787, 102 (2007). 

[28] A. N. Antonov et al, Nucl. lustrum. Meth. A637, 60 

(2011) . 

[29] T. Suda, "A construction proposal of an electron scat- 
tering facility for structure studies of short-lived nuclei" , 
Proposal for Nuclear Physics Experiments at RIBF 
NP1006 - SCRITOl (2010). 

[30] T. Suda, T. Adachi, T. Amagai, A. Enokizono, 
M. Hara, T. Hori, S. Ichikawa, K. Kurita, T. Miyamoto, 
R. Ogawara, T. Ohnishi, Y. Shimakura, T. Tamae, 
M. Togasaki, M. Wakasugi, S. Wang, and K. Yanagi, 
2012, 03C008 (2012). 

[31] E. Garrido and E. Moya de Guerra, Nucl. Phys. A650, 
387 (1999), arXivmucl-th/9902075 [nucl-th]. 

[32] E. Garrido and E. Moya de Guerra, Phys. Lett. B488, 
68 (2000), arXiv:nucl-th/0007024 [nucl-th]. 

[33] S. Ershov, B. Danilin, and J. Vaagen, Phys. Rev. C 72, 
044606 (2005). 

[34] Z. Wang and Z. Ren, Phys. Rev. C 70, 034303 (2004). 
[35] A. N. Antonov, M. K. Gaidarov, D. N. Kadrev, P. E. 

Hodgson, and E. Moya de Guerra, Int. J. Mod. Phys. 

E 13, 759 (2004). 
[36] A. N. Antonov, D. N. Kadrev, M. K. Gaidarov, 

E. Moya de Guerra, P. Sarriguren, J. M. Udi'as, V. K. 

Lukyanov, E. V. Zemlyanaya, and G. Z. Krumova, 

Phys. Rev. C 72, 044307 (2005). 
[37] C. A. Bertulani, Phys. Rev. C 75, 024606 (2007). 
[38] E. Khan, M. Grasso, J. Margueron, and N. Van Giai, 

Nucl. Phys. A800, 37 (2008), arXiv:0707.0363 [nucl-th]. 
[39] M. Grasso, L. Gaudefroy E. Khan, T. Niksic, 

J. Piekarewicz, O. Sorlin, N. Van Giai, and D. Vretenar, 

Phys. Rev. C 79, 034318 (2009). 
[40] Y. Chu, Z. Ren, T. Dong, and Z. W. Wang, Phys. Rev. 

C 79, 044313 (2009). 
[41] X. Roca-Maza, M. Centelles, F. Salvat, and X. Vifias, 

Phys. Rev. C 78, 044332 (2008), arXiv:0808.1252 [nucl- 
th]. 

[42] X. Roca-Maza, M. Centelles, F. Salvat, and X. Vinas, 
Phys. Rev. C 87, 014304 (2013), arXiv: 1207. 1679 [nucl- 
th]. 

[43] T. Dong, Y. Chu, Z. Ren, and Z. Wang, Phys. Rev. C 

79, 014317 (2009). 
[44] T. Dong, Z. Ren, and Z. Wang, Phys. Rev. C 77, 064302 

(2008). 

[45] J. Liu, Z. Ren, and T. Dong, Nucl. Phys. A888, 45 

(2012) . 



11 



[46] T. Dong, Y. Cliu, and Z. Ren, J. Phys. Conf. Ser. 381, 
012135 (2012). 

[47] C. Giusti, A. Meucci, F. D. Pacati, G. Co', and 
V. De Donno, Phys. Rev. C 84, 024615 (2011), 
arXiv:1105.1295 [nucl-th]. 

[48] G. Co', V. De Donno, P. Finelli, M. Grasso, M. An- 
guiano, A. M. Lallena, C. Giusti, A. Meucci, and 

F. D. Pacati, Phys. Rev. C 85, 024322 (2012), 
arXiv:1202.2512 [nucl-th]. 

[49] T. W. Donnelly, J. Dubach, and I. Sick, Nuclear Physics 

A 503, 589 (1989). 
[50] T. W. Donnelly and R. D. Peccei, Physics Reports 50, 

1 (1979). 

[51] C. J. Horowitz, Phys. Rev. C 57, 3430 (1998). 

[52] S. Abrahamyan et al. (PREX Collaboration), Phys. 
Rev. Lett. 108, 112502 (2012). 

[53] C. J. Horowitz, Z. Ahmed, C. M. Jen, A. Rakhman, 
P. A. Souder, M. M. Dalton, N. Liyanage, K. D. 
Paschke, K. Saenboonruang, R. Silwal, G. B. Franklin, 
M. Friend, B. Quinn, K. S. Kumar, D. McNulty, L. Mer- 
cado, S. Riordan, J. Wexler, R. W. Michaels, and G. M. 
Urciuoh, Phys. Rev. C 85, 032501 (2012). 

[54] PREX-II, Proposal to Jefferson Lab PAC 38, http:// 
hallaweb . jlab . org/parity/prex/prexll .pdf . 

[55] "CREX: parity-violating measurements of the weak 
charge distribution of ^**Ca to 0.03 fm accuracy". Pro- 
posal to Jefferson Lab PAC 39, http : //hallaweb . j lab . 
org/parity/prex. 

[56] O. Moreno, P. Sarriguren, E. Moya de Guerra, J. M. 
Udfas, T. W. Donnelly, and I. Sick, Nuclear Physics A 
828, 306 (2009). 

[57] S. Ban, C. J. Horowitz, and R. Michaels, Journal of 
Physics G: Nuclear and Particle Physics 39, 015104 
(2012). 

[58] X. Roca-Maza, M. Centelles, X. Vifias, and M. Warda, 

Phys. Rev. Lett. 106, 252501 (2011). 
[59] E. Moya de Guerra, O. Moreno, P. Sarriguren, and 

M. Ramon, Journal of Physics: Conference Series 366, 

012011 (2012). 

[60] D. Vretenar, P. Finelli, A. Ventura, G. A. Lalazissis, 
and P. Ring, Phys. Rev. C 61, 064307 (2000). 

[61] C. J. Horowitz, S. J. Pollock, P. A. Souder, and 
R. Michaels, Phys. Rev. C 63, 025501 (2001). 

[62] P. Finelh, N. Kaiser, D. Vretenar, and W. Weise, Nu- 
clear Physics A 735, 449 (2004). 

[63] P. Finelli, N. Kaiser, W. Weise, and D. Vretenar, 
"Nuclear density functional constrained by low-energy 
QCD," in Theoretical Nuclear Physics In Italy, edited 
by S. BofS, A. Covello, M. Di Toro, A. Fabrocini, 

G. Pisent, and S. Rosati (World Scientific, 2004) 
Chap. 5, pp. 45-52, arXiv:nucl-th/0411122 [nucl-th]. 

[64] P. Finelh, N. Kaiser, D. Vretenar, and W. Weise, Nucl. 

Phys. A770, 1 (2006). 
[65] C. Fuchs, H. Lenske, and H. H. Wolter, Phys. Rev. C 

52, 3043 (1995). 
[66] S. Typel and H. H. Woffer, Nuclear Physics A 656, 331 

(1999). 

[67] T. Niksic, D. Vretenar, P. Finelh, and P. Ring, Phys. 

Rev. C 66, 024306 (2002). 
[68] G. A. Lalazissis, T. Niksic, D. Vretenar, and P. Ring, 

Phys. Rev. C 71, 024312 (2005). 
[69] D. Vretenar, A. V. Afanasjev, G. A. Lalazissis, and 

P. Ring, Physics Reports 409, 101 (2005). 
[70] A. Meucci, F. Capuzzi, C. Giusti, and F. D. Pacati, 



Phys. Rev. C 67, 054601 (2003), arXivmucLth/0301084 
[nucl-th] . 

[71] A. Meucci, C. Giusti, and F. D. Pacati, Nuclear Physics 
A 739, 277 (2004), arXivmucl-th/0311081 [nucl-th]. 

[72] A. Meucci, C. Giusti, and F. D. Pacati, Nuclear Physics 
A 756, 359 (2005), arXivmucl-th/0501046 [nucl-th]. 

[73] A. Meucci, J. A. Caballero, C. Giusti, F. D. Pacati, 
and J. M. Udi'as, Phys. Rev. C 80, 024605 (2009), 
arXiv:0906.2645 [nucl-th]. 

[74] A. Meucci, J. A. Caballero, C. Giusti, and J. M. Udi'as, 
Phys. Rev. C 83, 064614 (2011), arXiv:1103.0636 [nucl- 
th]. 

[75] A. Meucci and C. Giusti, Phys. Rev. D 85, 093002 

(2012), arXiv: 1202.4312 [nucl-th]. 
[76] A. Meucci, C. Giusti, and F. D. Pacati, Phys. Rev. D 

84, 113003 (2011), arXiv:1110.3928 [nucl-th]. 
[77] A. Meucci, M. B. Barbaro, J. A. Caballero, C. Giusti, 

and J. M. Udi'as, Phys. Rev. Lett. 107, 172501 (2011), 

arXiv:1107.5145 [nucl-th]. 
[78] C. Giusti, A. Meucci, F. D. Pacati, G. Co', and 

V. De Donno, Journal of Physics: Conference Series 

366, 012019 (2012). 
[79] P. MoUer and J. R. Nix, Nucl. Phys. A536, 20 (1992). 
[80] A. Schiller, T. Baumann, J. Dietrich, S. Kaiser, W. Pe- 
ters, and M. Thoennessen, Phys. Rev. C 72, 037601 

(2005). 

[81] A. Meucci, C. Giusti, and F. D. Pacati, Phys. Rev. C 
64, 014604 (2001), arXiv:nucl-th/0101034 [nucl-th]. 

[82] A. Meucci, C. Giusti, and F. D. Pacati, Phys. Rev. C 
64, 064615 (2001), arXiv:nucl-th/0107052 [nucl-th]. 

[83] A. Meucci, Phys. Rev. C 65, 044601 (2002), arXiv:nucl- 
th/0111077 [nucl-th]. 

[84] M. Radici, A. Meucci, and W. H. Dickhoff, Eur. Phys. 
J. A 17, 65 (2003). 

[85] T. Tamae, Y. Sato, T. Yokokawa, Y. Asano, M. Kawa- 
bata, O. Konno, I. Nakagawa, I. Nishikawa, K. Hi- 
rota, H. Yamazaki, R. Kimura, H. Miyase, H. Tsubota, 
C. Giusti, and A. Meucci, Phys. Rev. C 80, 064601 
(2009). 

[86] C. Maieron, M. C. Martinez, J. A. Caballero, and J. M. 

Udi'as, Phys. Rev. C 68, 048501 (2003), arXiv:nucl- 

th/0303075 [nucl-th]. 
[87] A. Meucci, C. Giusti, and F. D. Pacati, Nuclear Physics 

A 744, 307 (2004), arXiv:nucl-th/0405004 [nucl-th]. 
[88] A. Meucci, C. Giusti, and F. D. Pacati, Acta Phys. 

Polon. B 37, 2279 (2006). 
[89] A. Meucci, C. Giusti, and F. D. Pacati, Nuclear Physics 

A 773, 250 (2006), arXiv:nucl-th/0601052 [nucl-th]. 
[90] J. A. Caballero, Phys. Rev. C 74, 015502 (2006), 

arXiv:nucl-th/0604020 [nucl-th]. 
[91] A. Meucci, C. Giusti, and F. D. Pacati, Phys. Rev. C 

77, 034606 (2008). 
[92] C. Giusti, A. Meucci, and F. D. Pacati, Acta Phys. 

Polon. B 40, 2579 (2009), arXiv:0906.2873 [nucl-th]. 
[93] J. A. Caballero, M. C. Martinez, J. L. Herraiz, 

and J. M. Udi'as, Phys. Lett. B 688, 250 (2010), 

arXiv:0912.4356 [nucl-th]. 
[94] A. V. Butkevich, Phys. Rev. C 82, 055501 (2010), 

arXiv: 1006. 1595 [nucl-th]. 
[95] A. V. Butkevich, Phys. Rev. C 85, 065501 (2012). 
[96] F. Capuzzi, C. Giusti, and F. D. Pacati, Nuclear 

Physics A 524, 681 (1991). 
[97] F. Capuzzi, C. Giusti, F. D. Pacati, and D. N. Kadrev, 

Annals of Physics (N.Y.) 317, 492 (2005), arXiv:nucl- 




0.05 - 



Figure 1. Neutron, panel (a), and proton, panel (b), distri- 
butions for the various oxygen isotopes we have considered. 
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Figure 3. Differential cross section for elastic electron scat- 
tering on ^^O at an electron energy e = 374.5 MeV and ^"Ca 
at £ = 496.8 MeV as a function of the scattering angle 9. 
Experimental data from [109] (^'^O) and [110] (^°Ca). 
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Figure 4. Differential cross section for elastic electron scat- 
tering on 14-28 Q ^ ^ 5 y^^y ^ function of 6. 
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Figure 5. Differential cross section for elastic electron scat- 
tering on ^"^"^^Ca at £ = 496.8 MeV as a function of 6. 
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Figure 6. Differential cross section of the reactions ^''0(e,e'), 
panel (a), and '"^Ca(e,e'), panel (b), for different beam ener- 
gies and electron scattering angles, e — 1080 MeV and = 32° 
for ^^0(e,e') and e = 841 MeV and 9 = 45.5° for *°Ca(e,e'), 
as a function of the energy transfer u. The RGF results are 
compared with the experimental data from [102] (^*0(e,e')) 
and [103] {*°Ca{e,e')). 
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CO [MeV] 

Figure 7. In panel (a) the differential RPWIA cross section for 
the inclusive QE (e,e') reaction on i^^^Sq a,t e = 1080 MeV 
and 6 = 32° is shown as a function of oj. In panel (b) the 
separate contributions of protons and neutrons are displayed. 
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Figure 8. Differential RGF cross section for the inclusive QE 
(e, e') reaction on ^''^^Sq ^j^g same kinematics as in Fig. 7. 
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Figure 9. Differential RPWIA cross section for the inclusive 
QE (e,e') reaction on ^'^"^^Ca at e = 560 MeV and 6 = 60° 
as a function of u}. 
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Figure 10. The same as in Fig. 9, but in the RGF model. 




Figure 11. Parity violating asymmetry parameter Ap„ for 
elastic electron scattering at e = 850 MeV as function of the 
scattering angle on ^''"'^^O. The green arrow emphasizes the 
evolution of A^-u as a function of the neutron number A*'. 
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Figure 13. First minimum positions of the asymmetry pa- 
rameter Apy as functions of A = (A'^ — Z)/Z for 1'*~28q 
36-56q^^ The black dashed hues represent the best linear fit. 




Figure 14. Panel (a): neutron density distributions lor some 
selected nuclei (^'^O and ^^O) that could be interpreted as 
candidates for "bubble" structure [39]. Panel (b): parity vi- 
olating asymmetry parameter Apv lor elastic electron scat- 
tering at £ = 850 MeV as a function of the scattering angle 

e. 




Figure 15. Panel (a): Theoretical weak charge density (red 
line) in comparison with the experimental error band as de- 
termined in Ref. [53] for ^"^Pb with the kinematics of the 
PREX experiment. The corresponding proton density profile 
is plotted with a black line. Panel (b): Asymmetry param- 
eter for ■^"^Pb as a function of the angle 6 (red line). The 
shaded area represents the interval covered by the acceptance 
function e{0) (see Ref. [52] for more details). The asymme- 
try parameter averaged over the acceptance (Apv) is in quite 
good agreement with the empirical value 0.656 ± 0.060(stat) 
± 0.014(syst) ppm. 
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Figure 16. Parity violating asymmetry parameters Ap^ for 
elastic electron scattering at e = 2.2 GeV as a function of the 
scattering angle 9 for ''^Ca, see Ref [55]. 




Figure 17. Parity-violating asymmetry for the quasielastic 
(e, e') reaction on O isotopes with the relativistic Green's 
function for e = 1080 MeV and 6* = 32°. 



